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with  a  particular  scaling  of  the  diffusion  coefficient  on  the  FEL  parameters.  Numerical  results, 
obtained  for  constant  radiation  amplitudes,  show  good  agreement  with  the  theoretical  predictions. 
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-highest  diffusion  rates  and  the  broad  continuous  spectrum  causes  the  lowest  diffusion  rates  under  con¬ 
stant  total  sideband  power.  The  diffusion  length,  measured  in  wiggler  periods,  is  independent  of  the 
beam  energy  y.  Serious  deterioration  in  the  FEL  efficiency  will  result  when  the  diffusion  length 
becomes  shorter  than  the  wiggler  length. 
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STOCHASTIC  ELECTRON  DETRAPPING  IN  FELS 
CAUSED  BY  SIDEBANDS 


I.  INTRODUCTION 


Multifrequency  effects  in  free  electron  lasers  (FEL)  become 
increasingly  important  as  progress  is  made  towards  high  power  operation. 
Growth  of  parasitic  frequencies  (sidebands^  "’)  has  been  predicted 
theoretically  and  has  been  observed  in  experiments^’^  as  well  as  in 
simulations^  ^  with  either  constant  or  tapered  wigglers^.  The 
efficiency  for  the  carrier  signal  is  reduced  and  the  optical  quality  is 
degraded  as  power  is  channeled  into  frequencies  apart  from  the  intented 
operation  frequency.  Another  potential  hazard  that  has  attracted  little 
attention  so  far  is  the  onset  of  chaotic  electron  motion  caused  even  by 
the  presence  of  a  single  frequency  sideband.  This  may  lead  to  extensive 
particle  detrapping  and  premature  loss  of  the  amplification  for  all  the 
radiation  modes  independent  of  frequency. 

So  far  the  attention  has  been  directed  to  the  linear  stability 

issue.  The  gain  for  small  sideband  signal  has  been  computed 

analytically^  ^  invoking  either  ensemble  averaging  over  single  particle 

trajectories  or  solutions  of  the  perturbed  kinetic  equation  for  the 

distribution  function.  Initial  results,  obtained  for  particles  localized 

near  the  bottom  of  the  ponderomot i ve  well,  and,  in  particular,  more 

2 

recent  results  including  all  trapped  and  untrapped  particles  with 
arbitrary  distributions,  have  made  it  clear  that  every  nontrivial 
distribution  dfQ/dJ  t  0  is  unstable  to  sideband  growth. 

Given  that  sidebands  cannot  be  eliminated,  the  growth  of  the 
unstable  modes  to  a  finite  amplitude  may  have  serious  effects  on  the 
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unperturbed  trajectories.  It  has  been  known  that  stochastic  behavior^ 

12 

is  an  intrinsic  property  of  perturbed  Hamiltonian  systems  .  Accordingly 
the  electron  motion  in  an  FEL  will  become  chaotic  when  the  sideband 
amplitude  exceeds  a  certain  threshold.  This,  in  turn,  will  result  into 
electron  detrappmg  and  potential  loss  of  amplification  for  the 
electromagnetic  fields. 

In  the  present  work  we  investigate  the  nonlinear  effects  caused 

by  sidebands.  The  threshold  for  stochastici ty,  above  which  unbound 

chaotic  motion  occurs  is  determined.  Once  the  stochastic  transition 

takes  place,  the  action  J,  a  constant  of  motion  in  the  unperurbed 

2 

system,  changes  in  a  random  manner.  The  ensamble  average  <AX  >  of  any 

physical  quantity  X  is  described  by  a  diffusion  equation.  Diffusion  of 

the  action  invariant  provides  a  measure  of  the  leakage  rate  across  the 

separatrix.  If  D  is  the  effective  diffusion  coefficient  in  action  space, 

2 

then  the  diffusion  length  =  Jg  /D  signifies  the  length  over  which 

approximately  half  of  the  deeply  trapped  particles  get  detrapped,  where 
is  the  action  at  the  separatrix.  We  show  that  a  single  frequency 
sideband  at  one  tenth  (1/10)  or  less  of  the  carrier  amplitude  suffices 
to  spread  irregular  motion  over  most  of  the  trapped  particle  domain. 
However,  given  that  the  interaction  time  of  an  electron  in  a  FEL  is 
short,  we  are  mainly  concerned  on  how  fast  this  diffusion  occurs.  The 
diffusion  rate  increases  and  the  diffusion  length  decreases  with 
increasing  sideband  ampli tude(s) .  Thus  a  critical  sideband  level  ac  can 
be  defined  above  which  the  diffusion  length  becomes  shorter  than  the 
wiggler  length  L  .  Obviously  the  power  level  for  the  sidebands  in  a  FEL 


must  not  exceed  ac,  otherwise  extensive  diffusion  and  premature 

detrapping  will  occur.  On  the  other  hand  enough  electrons  can  remain 

trapped  during  the  interaction  period  even  though  their  motion  has 

turned  stochastic  when  the  critical  amplitude  ac  is  much  larger  than  the 

threshold  for  stochasticity  a  . 

s 

A  clear-cut  relation  between  the  diffusion  rate  under  constant 
total  sideband  power  and  the  type  of  the  excited  sideband  spectrum  is 
discovered.  More  specifically  we  observe  three  regimes  in  the  simulation 
parameters  defining  the  sideband  spectrum,  corresponding  to  a  narrow,  a 
wide  discrete  and  a  wide  continuous  spectrum.  The  transition  from  one 
spectral  type  to  another  is  accompanied  by  an  abrupt  change  in  the 
diffusion  rates.  In  all  cases  we  find  the  diffusion  coefficient 
proportional  to  the  ratio  of  the  total  power  in  the  sidebands  to  the  FEL 
carrier  power.  The  coefficients  of  this  proportionality  depend  on  the 
spectral  type.  A  general  conclusion  is  that  the  diffusion  rate,  under 
constant  sideband  over  carrier  power  ratio,  decreases  with  increasing 
number  of  spectral  components.  The  diffusion  rate  for  a  single  sideband 
frequency  exceeds  that  of  a  broad  continuous  spectrum  by  orders  of 
magnitude  while  a  broad  discrete  spectrum  causes  intermediate  diffusion 
rates. 

The  analysis  also  shows  that  the  normalized  diffusion  coefficient 
does  not  depend  on  the  beam  energy  y  .  The  numerical  results  agree  well 
with  the  theory. 


II.  TRANSITION  TO  STOCHASTICITY 


We  consider  relativistic  electrons  streaming  along  the  z- 

direction  through  the  periodic  magnetic  wiggler  and  the  radiation  fields 

of  the  carrier  and  the  sideband.  We  assume  a  circularly  polarized  static 

wiggler  and  constant  amplitude  electromagnetic  waves.  We  assume  that  all 

waves  propagate  with  the  speed  of  light  c,  ignoring  the  small  correction 
2  2 

of  order  w  /w  from  the  dielectric  contribution  of  the  beam. 
P  r 

Electrostatic  contributions  to  the  fields  are  neglected  for  operation  in 

the  Compton  regime.  Both  the  wavenumber  k  (z)  and  the  amplitude  A  (z)  of 

w  w 

the  wiggler  may  change  slowly  in  z  on  a  scale  length  much  longer  than 
the  wiggler  wavelength  2 ri/k^, .  The  FEL  main  signal  wavenumber  kf  is 
doubly  Doppler  upshifted  from  the  wiggler  wavenumber  ky, 


k  =  2y  k  , 
r  z  w 


(1) 


2  -1/2 

with  the  upshifting  factor  yz  =  (1  -  0  )  and  =  «r/c(kr+kw) . 

-1  -1 

We  normalize  the  time  t  to  »  ,  the  length  z  to  k^  *,  the  mass 

2 

to  mg  and  the  vector  potentials  according  to  a^=|e|Aj/mec  where  the 
subscripts  w,  r,  and  s  stand  for  wiggler,  carrier  and  sideband 
respectively.  The  dimensionless  Hamiltonian  describing  the  electron 
motion  for  small  deviations  y  =  y  -  y  from  the  resonant  energy  y^  is 
given  by, 


WBWUfflHttWVX  WW  gTOtPTOOTTOg PgWCTrtWTOTWTO  WW  SKWDPTWnonPT  KKWjnam  tv  ni  rj  tv  r.’vvywv  vxviw.vw 


where  »J/r  and  5g  are  defined  by, 


d 

d^r 


k  a  a 
r  w  r 


sin^r  , 


(3) 


The  electron  motion  is  parametrized  by  the  traveled  length  z  inside  the 
wiggler.  The  term  sing/^  parametrizes  the  rate  of  change  for  the  resonant 
energy  caused  by  the  change  in  the  wiggler  wavelength  where  ^=0 
corresponds  to  an  untapered  wiggler.  The  term  Sg  in  the  sideband  phase 
is  the  Doppler  downshifted  difference  between  the  signal  and  the 
sideband  wavenumbers. 

In  the  absence  of  sidebands,  a  =0,  the  Hamiltonian  H  is 
integrable.  The  unperturbed  trajectories  in  the  ponderomotive  well  are 
given  by  H  (y^)  =  K  where  K  is  the  reduced  energy  in  the  ponderomotive 
frame.  These  orbits  take  the  simplest  possible  form  expressed  in  terms 
of  the  action  arg1?  variables  (J,  9),  defined  as; 


J  =  y(K,«I>),  9  =  fjJ  dg/  y(K,4/  )  ,  (4) 

where  K  =  Hq(J)  and  the  path  of  integration  is  over  the  unperturbed 
orbits.  Hamiltonian  (2)  is  transformed  under  the  canonical 
transformation  defined  by  Eq.  (4)  into, 


H(J,0;z)  =  Hq(J)  + 


a  a 
w  s 


J]  Q+(J)  cos(n9  +  6  z)  +  Q  (J)  cos(n9  -  &  z), 
Yr  n=0  n  s  n  s 

(5) 


Qn+(J)  are  the  Fourier  coefficients  obtained  by  the  decomposition  of  the 


perturbing  sideband  phase  <J/(J,0)  -  &sz  into  harmonics  of  the  angle  9, 
where  <p(J,  0)  is  obtained  by  inverting  Eq.  (4).  In  case  of  constant 
parameter  wiggler  J,  0  and  Qn  (J)  are  expressed  in  closed  forms  given  in 
Ref.  11. 

Hq(J)  is  independent  of  0  so  the  unperturbed  orbits  in  (J,  9) 
space  are  straight  lines  J  =  constant,  0  =  0Q  +  K^(J)z.  The  synchrotron 
wavenumber  K^J)  is  connected  to  the  bounce  length  and  the 
synchrotron  frequency  in  the  laboratory  frame  co^CJ)  with  the  relation 


Kb(J)  = 


dH  (J) 
o 


VJ) 


Lb<J) 


According  to  Eq.  (5)  new  resonances  emerge  when  a  sideband  is 
turned  on.  Defining  the  phase  of  the  nth  sideband  induced  harmonic 
=  n0  +  &sz,  the  stationary  phase  condition  reads, 


±  nK^(J)  -  &s  =  0,  or  ±n0z(A^(J)  -  8g  =  0. 


Thus  particles,  originally  in  unperturbed  orbits  J  =  J  ,  resonate  with 
the  sideband  when  the  nth  harmonic  of  their  synchrotron  period 
matches  the  downshifted  frequency  difference  between  the  sideband  and 
the  carrier  signal.  When  the  sideband  amplitude  becomes  finite  these 
resonant  particles  get  trapped  inside  secondary  ponderomotive  wells 
corresponding  to  the  different  harmonics  in  Hamiltonian  (4).  Thus  a 
single  frequency  sideband  causes  chains  of  secondary  islands  to  appear 


& 

I 

% 


I® 

$ 

i 


■iwV-Vvv 


inside  the  original  primary  island.  The  structure  of  the  phase  space  is 
shown  in  Figs,  la  -  Id.  They  are  surfaces  of  section,  created  by 
numerically  integrating  the  original  equations  of  motion  from 
Hamiltonian  Eq.  (2)  and  then  recording  the  intersection  point  of  each 
trajectory  with  the  plane  z  =  2 ti/ 6s .  The  y,  coordinates  are  plotted  on 
the  left  side  in  Figs,  la-ld.  The  plots  on  the  right  side  show  the  same 
surfaces  of  section  in  action-angle  variables,  produced  by  the 
transformations  (4). 

The  presence  of  just  one  sideband  frequency  suffices  to  transform 

the  regular  coherent  motion,  such  as  in  Figs,  la  and  lc,  to  the 

irregular  unbounded  motion  shown  in  Figs,  lb  and  Id  when  the  sideband 

amplitude  exceeds  a  certain  amplitude  ag  regarded  as  the  stochasticity 

threshold.  The  mechanism  for  this  radical  change  in  behavior  can  be 

briefly  described  as  follows.  The  trajectories  emanating  from  the 

unstable  fixed  points  (X-points)  of  a  secondary  island  corresponding  to 

a  given  harmonic  n  do  not  actually  join  smoothly  around  that  island. 

11  12 

They  intersect  infinite  times  with  each  other  ’  due  to  the  effect  of 

the  other  harmonics  n'  t  n,  creating  a  thin  layer  of  fussy  motion  that 

surrounds  each  island  chain.  As  the  amplitude  ag  increases  the  width  of 

each  island  increases  and  so  does  the  thickness  of  the  stochastic  layer 

around  that  island.  At  a  given  point  the  stochastic  layers  around  the 

1 3 

two  neighboring  island  chains  n  and  n+1  overlap  ,  allowing  particles  to 
hop  from  one  island  to  another.  This  signifies  the  beginning  of 
unbounded,  random  motion  in  J  characterized  as  stochastic  diffusion. 


v'~-w-sV- 


The  overlapping  condition 


that  defines  the  transition  to 
11  13 

stochastici ty  is  usually  expressed  by  ’ 

SJ  +  SJ  ,  >  \  AJ  ,  (8) 

n  n+1  3  n  v 

where  SJ  ,  SJ  .  are  the  separatrix  half-widths  and  AJ  =  J  ,  -  J  is 
n  n+1  r  n  n+1  n 

the  distance  between  the  separatrix  centers  for  the  n  and  n+1  harmonics 

respectively.  SJn  and  AJn  are  evaluated  from  Hamiltonian  (5)  using  the 

familiar  methods  of  nonlinear  dynamics  in  terms  of  the  sideband 

amplitude  a^  and  the  shear  dK^(J)/dJ.  The  detailed  computation, 

performed  in  Ref.  14,  shows  that  the  threshold  for  stochastici ty  a 

depends  primarily  on  the  frequency  separation  (oo^  -  and  less  on 

the  other  FEL  parameters.  The  ratio  a  /a  varies  in  a  semioscillatory 

s  r 

manner  between  typical  values  of  0.1  and  1.  The  above  threshold  was 
derived  for  a  single  sideband  frequency.  Mult i frequency  sideband  spectra 
generally  tend  to  reduce  a^. 

III.  NARROW  SPECTRUM  DIFFUSION 

We  examine  first  the  diffusion  caused  by  a  narrow  sideband 
spectrum,  meaning  that  the  width  of  the  sideband  spectral  line  is  much 
less  than  its  distance  from  the  carrier  frequency.  We  model  this  case  by 
a  single  frequency  sideband.  Figures  lb  and  Id  show  a  typical  phase 
portrait  when  the  sideband  amplitude  a^  is  slightly  above  .  Two 
different  kinds  of  regions  coexist:  a  stochastic  Legime  where  diffusive 


behavior  prevails,  interrupted  here  and  there  by  islands  of  regular 
motion,  remnants  of  the  original  regular  motion.  The  stochastic  regimes 
are  interconnected  allowing  particle  transport,  hence  ag  was  defined  as 
global  stochastici ty  threshold. 

Figure  2  illustrates  the  stochastic  transition  by  plotting 
selected  particle  trajectories  J(z)  with  z  measured  in  wiggler  periods. 
The  sideband  amplitude  increases  from  well  below  in  Fig.  2a  to  well 
above  in  Fig.  2c.  The  initial  conditions  for  the  particles  remain  the 
same  in  all  cases.  The  oscillations  in  Fig.  2a  have  the  period  of  a 
particle  trapped  around  a  secondary  island. 

When  ag  is  increased  well  above  ag  the  chaotic  motion  engulfs 
almost  100%  of  the  phase  space.  Total  stochastization  of  the  island 
interior  occurs  roughly  when  the  sideband  amplitude  grows  to  the  point 
where  the  stable  fixed  point  y  =  0,  =  n  +  ^  at  the  centre  of  the 

original  island  becomes  unstable. 

The  two  questions  of  practical  interest  are:  (a)  what  percentage 
of  the  particles  will  eventually  get  detrapped  and  (b)  how  fast  do  they 
leak  outside  the  separatrix.  For  our  uniform  initial  distribution  the 
maximum  fraction  of  particles  becoming  detrapped  equals  the  fraction  of 
the  inside  the  separatrix  area  that  becomes  chaotic.  In  Fig.  3a  we  plot 
the  fraction  f^  of  the  particles  that  cross  the  original  separatrix 
as  a  function  of  the  traveled  wiggler  length  for  values  of  a^a^.  below 
the  threshold  for  extended  stochasticity.  In  all  cases  an  initial  stage 
of  quick  diffusion  is  followed  by  a  long  period  where  the  average  number 
of  untrapped  particles  remains  practically  constant.  The  results  are 


consistent  with  the  existence  of  a  boundary  in  phase  space  (RAM  surface) 
separating  two  regimes:  the  one  of  unbounded,  chaotic  motion  from  the 
one  filled  with  regular,  coherent  orbits  of  particles  that  remain 
trapped.  Only  electrons  in  the  area  between  the  last  integrable  surface 
and  the  old  unperturbed  separatrix  will  diffuse  until  that  area  is 
depleted.  A  fraction  1  -  f^  of  the  original  primary  island  area  will 
remain  trapped  for  an  arbitrarily  long  time,  as  long  as  ag  remains  below 
the  threshold  associated  with  the  particular  sideband  frequency.  This 
fraction  is  shrinking  as  ag  increases  and  the  bucket  "peels  off".  The 
situation  when  a  exceeds  a  is  shown  in  Fig.  3b.  The  fraction  of 
untrapped  particles  f^  increases  monotonically ,  reaching  1  in  some 
cases,  meaning  complete  absence  of  particle  confinement  in  the  bucket. 
All  particles  can  eventually  escape  with  a  rate  that  increases  with 
increasing  as-  The  length  over  which  approximately  half  of  the  initially 
trapped  particles  get  detrapped  will  be  discussed  in  the  next  section, 
in  comparison  with  the  diffusion  rates  from  other  types  of  sideband 
spectra.  It  is  evident  from  the  results  in  Fig.  3  that  the  diffusion 
length  at  constant  amplitude  ag  depends  on  the  sideband  frequency. 

IV.  BROAD  FREQUENCY  BAND  DIFFUSION 

So  far  stochastic  electron  detrapping  caused  by  a  single 
frequency  sideband  has  been  examined.  Here,  we  consider  the  situation 
when  many  sideband  frequencies  have  been  excited.  We  will  make  a 
distinction  between  a  continuous  and  a  discrete  spectrum.  In  case  of  a 


■&A'M  »■>:>.»  1,1  li>  Q  l.ri 


discrete  spectrum  the  distance  between  two  nearby  sideband  frequencies 
is  much  larger  than  the  width  of  an  individual  spectral  line.  In  the 
opposite  case  when  various  peaks  in  the  spectrum  merge  together  we  will 
talk  about  a  continuous  spectrum.  We  may  model  numerically  both  cases  by 
introducing  a  modulation  in  the  sideband  phase  of  Hamiltonian  Eq.  (4), 


k  ~0  a  a  a  a 

w  z  w  r  w  s 

H(y,<J»;z)  =  —  y  +  - (cos<J/  +  <J/  sim|/  )  +  - cos(\p  +  A  sin  vz 

Yr  Yr  Yr 


that  is  transformed  in  action-angle  variables  as, 


&sz)’ 


H(J , 0; z)  =  Hq(J) 


a  a  <*> 

W  S  r 


E  Jm(A>  E  Qn<J)  cos[n0  +  Ss(m)z]  +  Qn(J)  cos[n0  -  6s(m)z]  . 

=-°°  n=0 


r  m=-°° 


The  frequency  mismatch  values  6  (m)  and  the  corresponding  sideband 

s 

frequencies  <*>s(m)  are  given  by, 


6  (m)  =  S  +  m  v,  w  (m)  =  «  +  2  m  y  v, 

sv  '  o  ’  sv  so  z 


where  w  =  w  +  2y  6  .  Since  the  Bessel  function  coefficients  become 
so  r  z  o 

vanishingly  small,  Jm(A)  <<  1  for  A  >  m,  the  width  of  the  spectrum  is 

2 

given  by  D&s  ~  A  v  or  D«s  =  2yz  Av. 

Ir.  order  to  examine  the  connection  between  diffusion  rates  and 
the  types  of  the  sideband  spectra,  we  divide  the  latter  into  three 
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general  categories:  narrow,  broad  discrete  and  broad  continuous.  The 
passage  from  one  regime  to  the  other  is  characterized  by  abrupt  changes 
in  the  diffusion  coefficients.  Both  cases  of  the  broad  spectrum  are 


characterized  by  a  width  D«  in 
1  s 

than  the  upshifted  synchrotron 
equivalently, 


the  excited  frequencies  that  is  larger 

2 

frequency  D (*>s  >  2y z  u^,  or 


w. 

v  >  f-  ,  (12) 

with  A  >>  1.  The  further  distinction  between  discrete  or  continuous 

spectrum  is  related  to  the  separation  between  nearby  frequencies.  We 

1/2 

find  that  when  w^/A  >  \>  >  m^/A  the  diffusion  rate  agrees  well  with 

the  quasilinear  diffusion  coefficient.  A  different  coefficient  is 

1/2 

derived  for  the  case  when  \>  >  w^/A  >  w^/A,  in  agreement  with  the 

1/2 

numerical  simulations.  Consequently,  the  separation  w^/A  between 
nearby  modes  marks  the  transition  from  a  discrete  to  a  continuous  type 
of  behavior. 

(a).  Broad  Discrete  Spectrum. 

We  now  evaluate  the  diffusion  coefficient  for  a  broad,  discrete 

1/2 

spectrum  associated  with  frequency  separation  v  >  w^/A  .  Due  to  the 
presence  of  many  frequencies  in  the  spectrum  J(z)  executes  a  complicated 
oscillatory  motion  with  the  average  <  J  >  changing  very  little  most  of 
the  time.  J  however  receives  a  large  kick  AJ  near  resonances,  where  the 


phase  $  =  t  +  A  sin  \>z  ±  §sz  of  the  multifrequency  Hamiltonian  (9) 

varies  slowly.  The  resonant  condition  is 

k  y 

+  A  v  sin  vz.  +6  =  0,  (13) 

dz  y  l  -  s 

r 

at  some  z  =  z^.  Given  that  kvy/Yr  <  collective  effects  due  to  many 

frequencies  are  important  for  the  resonance  in  Eq.  (13)  when  A  \>  >  «.q. 

On  this  basis  inequality  (12)  signifies  the  transition  from  a  narrow  to 

a  broad  spectrum.  Let  us  consider  the  case  A  v  >>  Then  the 

resonances  occur  at  z^  =  in/v,  i  integer,  and  the  interval  between 

successive  resonances  is  Az  =  n/v.  Expanding  the  phase  $(z)  in  the 

equation  of  motion  for  J  around  the  resonance  z^,  and  extending  the 

2  2 

limits  of  the  z  integration  from  to  +“  we  obtain,  for  A  \>  >>  , 


a  a  (  2n  11/2  , 

DJi  “  —2  («„  *  ?)  lV<*„x>  -  V<*i>  I  - 

Y  \  A  v  / 


where  V(t)  =  cost  +  t  sintr.  The  resonant  phases  between  two 
successive  jumps  AJ  become  quickly  decorrelated  when  ag  grows  above  the 
stochasticity  threshold.  Assuming  complete  decorrelation  in  between 
two  successive  jumps  we  obtain, 
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where  <....>  is  the  ensemble  average  over  For  practical  purposes 
it  is  more  convenient  to  rescale  the  diffusion  coefficient  so  that  the 
distance  z  =  z/X^  is  measured  in  terms  of  wiggler  wavelengths  and  the 
action  J  =  J/Js  signifies  the  location  relative  to  the  separatrix.  In 
these  units,  using  Eq.  (4)  for  Jg  and  setting  v  /  <*^(0)  =  r  we  obtain 
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w  J 
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U  ♦  aw 


(16) 


The  term  g  is  a  scaling  factor,  the  ratio  of  the  untapered  wiggler 

separatrix  to  that  of  a  tapered  wiggler,  depending  only  on  ^  ,  g  = 

Js('l'r=0)/,'Js^r^ '  T^e  lyP*03!  diffusion  length  L^,  the  traveled  distance 

inside  the  wiggler  over  which  the  average  trapped  particle  crosses  the 

2  2 

separatrix  is  estimated  from  <AJ  >  =  J  =  D  L, .  Thus,  the  diffusion 

s  w  d 

length  in  wiggler  periods  N,  =  L,/X,  is  the  inverse  of  D  ,  N,  =  1/D  . 

daw  — w  d  ~v 


(b).  Broad  Continuous  Spectrum 


Next  we  consider  the  case  of  a  sideband  wave  package, 
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One  condition  for  the  applicability  of  the  quasilinear  theory  is 


that  the  phase  mixing  due  to  Dkg  occurs  much  faster  than  the  bounce 


period  around  a  secondary  island  in  phase  space.  This  way  electrons, 


that  otherwise  would  execute  periodic  orbits  around  some  fixed  point, 


lose  coherency  fast  enough  and  a  random  motion  of  the  Fokker-Planck  type 


sets  in.  This  condition  is  satisfied  for  spectral  widths  Dkg  of  order, 


O  V  v  ^  V 

nr >?  2  Yz  — 

s  k 


2  Kb  (0) 
,  2 


The  overlapping  condition  among  nearby  island  chains,  similar  to  Eq. 


(8),  must  also  be  satisfied.  The  threshold  for  overlapping  is  much  lower 


than  ag  related  to  the  single  sideband  frequency.  Then  the  evolution  of 


the  initial  distribution  f q ( J )  is  globally  described  by  a  diffusion 


equation, 


M  .  3_  D  (J)  i f 

dz  ~  aj  qWJ  aj 


Applying  the  standard  procedures  of  the  quasilinear  theory  and  taking 


the  limit  of  small  growth  rate  for  the  sidebands,  Im(kg)/kg  <<  1,  we 


obtain  the  quasilinear  diffusion  coefficient, 


VJ>  = 


k  a2y2V(k  )k  J 
r  v  z  -s  so  w 


V“b'J) 


We  express  D^(J)  in  normalized  units,  with  the  wiggler  wavelength  as 
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the  unit  length  and  the  action  J  at  the  separatrix  as  the  unit  action. 

s 

Chosing  the  value  J  =  J  /2  inside  D  (J)  we  obtain  an  estimate  for  the 

s  q 

effective  diffusion  coefficient  associated  with  the  uniformly  filled 
distribution, 


D 

-q 


2n  VJ)  _  “  g  avar  Vs(t)  ? 

kw  Js2  4  A  (l+aw2)  Dkg  Vr(t) 


(20) 


where  g  is  the  same  scaling  factor  as  in  Eq.  (16). 

Note  that  both  expressions  (16)  and  (20)  for  D,  corresponding  to 
the  two  different  spectral  types,  are  independent  of  rr*  1'hus,  for  the 
same  wiggler  parameters  and  total  sideband  power,  the  detrapping 
distance  in  wiggler  periods  is  independent  of  the  electron  beam  energy. 

A  more  detailed  treatement  of  the  subject,  including  the 
derivation  of  the  formulas  (16)  and  (20)  is  found  in  Ref.  14. 


V.  NUMERICAL  RESULTS 


The  trajectories  J(z)  on  the  left  side  of  4a-4c  are  generated  by 

the  same  initial  conditions  for  the  electrons,  the  same  FEL  parameters 

2 

a  ,  a  and  k  ,  and  the  same  averaged  sideband  power  <a  >.  The  spectral 
w  r  w  °  s 

parameters  A  and  v  however  are  different  so  that  each  of  the  cases  (a) 
(b)  and  (c)  corresponds  to  one  of  the  three  spectral  types  defined 
earlier.  The  dashed  line  marks  the  position  of  the  unperturbed 
separatrix  Jg.  Without  the  sidebands  the  trajectories  would  be  straight 


lines.  The  corresponding  distribution  functions  f(J,z)  are  plotted  on 

the  right-hand  side  at  the  beginning,  z  =0,  halfway  inside,  z  =  50X  , 

and  at  the  end,  z  =  100X  of  the  wiggler. 

w 

The  numerically  computed  diffusion  length  in  wiggler  periods 

=  1/D  is  plotted  against  the  sideband  to  carrier  power  ratio  P  =  Z 
2  2 

as  ^wn^'/ar  =  ^s^r  *n  ^*8"  ^  f°r  t^ie  three  different  types  of  spectra  . 

We  have  integrated  numerically  the  equations  of  motion  for  400  particles 

of  a  uniform  initial  distribution  inside  the  bucket.  The  field 

intensities  remained  constant  at  a  =  5x10  a  =2  and  y  =  25.  A 

r  w  r 

clear  separation  in  the  diffusion  rates  is  observed  among  the  various 
spectral  types.  The  narrow  frequency  results  (triangles)  were  obtained 
using  the  Hamiltonian  (2)  with  a  single  sideband  frequency  wg/wr  = 
1.016.  The  results  for  a  broad  discrete  spectrum  (circles)  were  obtained 
using  (9)  with  A  =  20,  oog/wr  =  1.016  and  v  =  0.5  6g.  The  continuous 

spectrum  (squares)  was  modeled  by  A  =  100,  v  =  0.05  6g.  The  solid  lines 

correspond  to  the  theoretical  results  of  Eqs.  (16)  and  (20),  in  good 
agreement  with  the  numerics.  Theoretical  predictions  for  the  single 
frequency  case  were  not  made.  We  stress  however  the  difference  between 
single  frequency  results  and  quasilinear  theory  in  this  case.  The 
agreement  that  has  been  observed  in  some  other  cases  is  not  generic  but 
particular  to  certain  systems. 

Finally,  in  Fig.  6  we  plot  the  diffusion  coefficient  for  a 
uniformly  filled  bucket  as  a  function  of  the  energy  y^ ,  fixing  the 
wiggler  parameters.  It  is  clear  that  the  diffusion  rate  (measured  again 
in  number  of  wiggler  periods)  is  independent  of  the  beam  energy, 


provided  the  synchrotron  frequency  stays  in  the  same  parameter 

regime. 

Once  the  diffusion  coefficients  are  known  some  estimate  can  be 
made  of  the  related  reduction  in  efficiency  over  the  wiggler  length. 


Assuming  that  n^(z)  =  n^(0)  exp(-z/L^)  the  number  of  detrapped  particles 


between  z  and  z  +  Az  is  An(z)  =  n^(0)L^  exp(-z/L^)  A z.  These  particles 


gave  up  an  amount  of  energy  AE(z)  =  [ vr ( 0 )  -  Yr(z)]  An(z)  as  radiation. 


Integrating  AE(z)  over  the  wiggler  length  for  a  linearly  tapered  wiggler 
rr(z)  =  Yr(0)  -  Ay  z  /  Ly  we  find  the  total  energy  extracted  from  the 
particles  that  were  detrapped  at  some  point  inside  the  wiggler.  Adding 
the  contribution  [Yr(0)  -  Yr(Lv)]  n^(Ly)  from  the  particles  that 
remained  trapped  throughout  the  wiggler  length  we  come  up  with, 


"<L»>  ■  no  r  (  1  -  eKp  <-  lt>  )■ 

w  d 


(21) 


where  hQ  =  Ay/y,.(0)  is  the  efficiency  without  induced  diffusion.  The 


loss  of  amplification  will  in  general  be  distributed  among  all  the 
radiation  modes  and  (21)  reflects  the  total  power  loss  in  all 
frequencies . 
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Figure  1. 


Plots  of  surfaces  of  section  in  both  y»  'P  and  J,  0 

representations,  illustrating  transition  to  chaotic 

behavior.  The  parameters  are  ay  =  2,  a^  =  5x10 and 

initial  y^  =  25.  Figures  (a)  and  (b)  correspond  to  tapered 

wiggler,  <|>r  =  4ri/3,  with  «s/(*>r  =  1.016  and  (a)  ag  =  2x10” 

^ ,  (b)  ag  =  lxlO-^.  Figures  (c)  and  (d)  correspond  to 

untapered  wiggler  with  ws/“r  =  1-024  and  (c)  ag  =  2xl0-^, 

(d)  a  =  1x10  . 
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Parameters  are  a„  =  2,  a  =  5x10  ,  y  =25  and 
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Figure  4.  Particle  response  under  different  sideband  spectra.  On  the 

left-hand  side  we  plot  the  action  J  against  the  number  of 

wiggler  periods  N  for  selected  particles.  On  the  right- 
w 

hand  side  we  plot  the  corresponding  distribution  function 
f(J)  at  =  0,  50,  and  100  .  In  all  cases  the  total 
sideband  fractional  power  P  =  Vg/Wr  =  0.36  and  ay  =  2,  a^ 
=  5x10  ,  Yr  =  25.  (a)  corresponds  to  a  discrete  sideband 

spectrum  with  o^/c*^  =  1.016.  (b)  corresponds  to  a  wide 

discrete  spectrum  and  (c)  to  a  wide  continuous  spectrum. 
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1/D  as  a  function  of 


The  e-folding  diffusion  length  = 

the  sideband  fractional  power  P  =  W  /W  for  a  =  2,  a  = 

r  s  r  w  ’  r 

5x10  Yr  =  25  and  ws/wr  =  1-016.  Squares  correspond  to  a 
continuous,  dots  to  a  wide  discrete  and  triangles  to  a 
single  frequency  spectrum.  The  upper  and  lower  solid  lines 
correspond  to  the  theoretical  results  Eqs.  (16)  and  (20) 
respectively. 
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Plots  of  the  normalized  diffusion  coefficient  D  against 

_4 

the  electron  energy  y  for  a  =  2,  a  =  2x10  and  a  = 
7.5x10 Squares  correspond  to  a  continuous,  dots  to  a 
wide  discrete  and  triangles  to  a  single  frequency 
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